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Abstract 

It is believed that grain boundary segregation of light interstitials can serve as the possible mechanism of thermal 
stability in commercially pure nanostructured titanium alloys. In this paper, using first-principles calculations, we 
show that independent segregation of C, N and O atoms at 27 high angle grain boundary in a-T\ is energetically 
unfavourable. The presence of interstitial elements near the grain boundary plane results in the increase of the grain 
boundary width and specific formation energy. 
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1. Introduction 

An increased (decreased) concentration of impuri- 
ties at grain boundaries (GB), known as segregations 
can considerably influence various properties of poly- 
crystalline materials |T]. Besides causing an embrittle- 
ment effect, they can increase grain boundary cohesion 
and significantly improve thermal stability of nanocrys- 
talline materials 12]. Hence, it is important for the 
design and optimisation of nanomaterials, which show a 
large specific area of interfaces that can give rise to seg- 
regations, to understand the influence of segregations on 
materials properties. 

Commercially pure nanostructured a-titanium (CP- 
nTi) alloys with hep lattice structure have a high bio- 
compatibility and superior strength in comparison to 
coarse-grained Ti lllj 

la. 

Being a promising material 
for medical lifetime dental implants the problem of ther- 
mal stability is very important for CP-nTi. There are no 
data in literature about long term behaviour of this mate- 
rial, because the mechanisms that are thermally stabilis- 
ing nTi are not understood yet. Due to a technologically 
necessary heat treatment of implants during production, 
an improvement of their microstructure stability is a sig- 
nificant task, which needs first of all a detailed under- 
standing of existing mechanisms of grain boundary sta- 
bilisation. CP titanium alloys (Grade 1-4) contain a suf- 
ficiently large amount of light element impurities such 
as C, N and O Usll . Hence, two most probable mecha- 
nisms of thermal stability are possible in nTi: (i) grain 
boundary pinning by precipitates (dispersed particles of 
carbides, nitrides, etc. formed from solid solution) and 



(ii) lowering of driving force for grain growth and re- 
duction of grain boundary mobility due to segregation 
of existing impurities |2]. The possibility of formation 
of Ti-C precipitates in a-T\ was considered in our recent 
work OlSll . The purpose of the current work is to con- 
tinue the investigation of thermal stability mechanisms 
in CP-nTi, considering the case of light impurity seg- 
regation at a-Ti grain boundaries in comparison with 
previous results of Ti-C particles formation. 

Semenova et al. flB'] have recently observed an in- 
creased concentration of C, N and O at grain boundaries 
in CP-nTi by atom probe tomography, confirming that 
segregation of light elements is possible at a-Ti grain 
boundaries. However, one could not derive clear con- 
clusions about their influence on materials properties, 
since many factors are acting together Does the segre- 
gation correspond to an equilibrium state? Does a co- 
segregation of elements take place? Is the formation 
of clusters at grain boundaries favourable? How strong 
is the interaction of impurities with the GB? Answer- 
ing these questions solely with experimental methods 
is very hard or even impossible and requires the use of 
computer simulation [17]. Recently, studies of segrega- 
tion from the first principles have become more com- 
mon il8l433ll . however, in the case of hep metals such 
activity is almost absent. In this work we want to fiU 
this gap with the example of hep titanium. 

In summary, a systematic study of segregation in 
Ti will allows taking further steps towards discover- 
ing the most effective mechanism of thermal stability 
in CP-nTi. One may note that quantitative studies of 



thermal stability of grain boundary structures, includ- 
ing kinetic details of mechanisms are most efficiently 
done by molecular dynamics and Monte Carlo simula- 
tions. However, these methods are based on empirical 
interatomic potentials, and therefore a benchmarking 
with first principles calculations (which are currently 
not available) is also in this case decisive iHtIi . The basis 
for such calculations is the determination of segregation 
formation energies, which are needed for impurity atom 
to diffuse from a bulk site to a GB site. 

Thus, in this paper we investigate interaction of C, 
N and O impurities with E7[0001](1230) GB in a-Ti 
from first principles. The results include segregation 
energies of impurities for different positions near and 
at GB plane. The influence of impurities on the grain 
boundary is accounted by the relaxation of GB atomic 
structure. We involve analysis of the electronic structure 
to explain several features of interaction between grain 
boundary and impurities. 

2. Details of calculation 

Calculations of the full energies and optimized ge- 
ometries were performed in the framework of the 
density-functional theory (DFT) lf34l 35 1 within the 



Table 1: Specific excess energy y (J/m-) of grain boundaries before 
and after relaxation, sizes a, b and c (A) of considered supercells, and 
distance rfxi-Ti (A) between Ti atoms in [1010] direction. 



generalized gradient approximation (GGA) using the 
Perdew-Burke-Wang |36| functional and projected aug- 
mented wave (PAW) method l37] (ABINIT Q). 
We considered the following valence electronic states: 
3i, 3/7,4i, 3c/ for Ti and 2s,2p for C, N and O. To 
take off any restrictions during relaxation, we set num- 
ber of point group symmetry operations to one in all 
cases. The calculations were performed under three- 
dimensional periodic boundary conditions. An or- 
thorhombic supercell containing 56 Ti atoms was con- 
structed using calculated in our previous work [15] lat- 
tice constants of hep Ti. The volume of cell was sub- 
sequently relaxed to account the influence of cell sizes 
and refine the values of lattice constants. The commen- 
surate supercell containing 54 Ti atoms was eniployed 
as a starting model for the S7[0001](1230) GB jMEj 
(Fig. [T]). In order to find the optimal geometry of GB, 
the supercell was relaxed with fixed dimensions (as for 
bulk Ti) in the plane of grain boundary. 

The sizes of supercells and lengths c/xi-Ti between Ti 
atoms in [1010] direction (coincide with a lattice con- 
stant in the case of bulk Ti) are listed in Table [T] The 
value of c/xi-Ti for bulk Ti supercell is in agreement with 
experimental lattice constant. In the case of supercell 
with GB the variation of t/xi-ii is caused by the influ- 
ence of grain boundaries and related to the finite sizes of 
cell. However, recent studies confirm that selected cell 
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sizes are enough for the description of grain boundaries 
structure in hexagonal materials [41 [ and GB segrega- 
tion energies [33il . 

To study the interaction of impurities with the GB, 
we placed atoms at different positions in optimised su- 
percell and relax it according to all internal coordinates 
and size of supercell in the direction perpendicular to 
the grain boundary plane. The calculations were per- 
formed by using an energy cut-off of 540 eV for the 
plane-wave basis set. The Brillouin-zone (BZ) inte- 
grals were approximated using the special k-point sam- 
pling of Monkhorst and Pack [42] with 2x1x4 
grid. We have used the Methfessel-Paxton 143] smear- 
ing for Brillouin-zone integration with a smearing width 
of 0.027 eV. Such relatively small value does not in- 
fluence the results but ensures faster convergence. The 
structural optimization was performed until the forces 
acting on each atom became less than 25 meV/A. The 
computational setup ensures that differences in the seg- 
regation energies of interstitials and formation energies 
of Ti grain boundaries are converged to within 25 meV 
and 0.05 J/m^, respectively. The method of calculation 
and Ti-C PAW potentials were checked by computation 
of lattice and elastic constants of several Ti and Ti-C 
phases in our previous study [15.1 . 



3. Results and Discussion 

3.1. Structures and energies of the pure GB 

There are periodic and free boundary conditions 
(PBC and FBC) which are mainly used for grain bound- 
ary modelling. Free boundary conditions allow to ex- 
amine one grain boundary in simulation cell at the cost 
of adding two open surfaces. This is useful for the 
cells with two non-equivalent grain boundaries. How- 
ever, there is significant influence of open surfaces on 
atomic structure of supercell. The attempt of reduction 
of such influence by introduction of fixed layers cause 
additional restrictions on the relative moving of adjacent 
grains during relaxation. In the present work, despite 
the existence of two non-equivalent grain boundaries in 
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simulation cell, to avoid difficulties related to open sur- 
faces, we use periodic boundary conditions. This is rea- 
sonable due to the specific atomic configuration of cell 
with one fully coherent grain boundary obtained after 
relaxation. 

As initial structure for grain boundary the ideal sym- 
metric coincidence-site lattice (CSL) grain boundary 
S7[0001](1230) with ^csl=21.8° is used. This bound- 
ary perfectly agrees with experimentally observed GBs 
in a-Ti ll44ll . The procedure of construction of this GB is 
described elsewhere (3^. After the mathematical con- 
struction of supercell there are pairs of atoms that lie too 
close to each other. Hence, according to 1I39I1 each pair 
was replaced by one atom. 

Two simulation cells are shown in Fig. [TJa) to indi- 
cate the structure of both grain boundaries. Ti atoms 
in the bulk regions and at GBs are shown with different 
colours to guide the eye. The atomic configuration at Gl 
and G2 grain boundaries (see Fig. [IJa)) in both layers 
compose of pentagons, which is in consistent with min- 
imum binding energy structure obtained in [39|. How- 
ever, there are differences between structures of Gl and 
G2 in the bottom layer A (see Fig. [TJa))- The pentagon 
in this layer at G2 is stretched in [1230] direction rel- 
ative to pentagon at Gl grain boundary. Moreover, the 
upper pentagon at G2 is shifted relative to the bottom 
pentagon in [5410] direction compared to Gl structure. 
These differences arise from the fact that the coincide 
site lattice exists only in layer B for the given rotation 
angle. 

Two simulation cells after relaxation are shown in 
Fig. [TJb). The main difference from non-relaxed struc- 
ture is that grains are shifted relative to each other along 
the [0001] direction. The interesting feature is that shift 
occurs only at Gl plane and is absent at G2. The values 
of shifts Ari in [5410] direction and Ara in [0001] di- 
rection are averaged by several atoms in the bulk region 
of grains and calculated relative to the initial state. The 
shift is almost absent in [5410] direction for the pure 
relaxed GB, while the value of Ar2 is 0.45 A. 

The atomic structure of Gl and G2 grain boundaries 
was also significantly changed after relaxation. In the 
case of Gl grain boundary, T3 and T4 atoms are shifted 
in opposite directions breaking the symmetry of pen- 
tagon and forming more close-packed configuration that 
is more preferable for titanium. The reduce of cell size 
in [1230] direction due to the relaxation is 3.7 % (see. 
TablelB- 

The atomic structure of G2 is highly symmetrical. 
The visual analysis allows to determine w-Ti phase 
1 4511 at G2. The conventional cell of w-Ti carved out 
from G2 is shown in Fig. |2] where marked Ti atoms 



matches with Fig. [T] The orientation relationship is 
(OOOIX, II (1210)„ and [1230]„ || [0001]^, which is 
agreed with Silcock [46] relationship for planes, but 
differs for directions ([1120]^ || [0001]^ for Silcock). 
We did not find experimental confirmations of such 
[1230]q. II [0001](j relationships for directions, but the 
obtained cj-a interfaces are quite coherent. The lattice 
parameters of w-Ti phase at grain boundary are in agree- 
ment with that for bulk w-Ti calculated in Iil5il . The in- 
terface period in [0001] direction (4.62 A) matches very 
well with a parameter of bulk w-Ti (4.58 A). The period 
in [5410] direction (7.77 A) matches well with V3a in 
bulk w-Ti(7.92 A). 

We study segregation at Gl grain boundary using the 
PBC, because the influence of the coherent G2 grain 
boundary is not larger than that of the open surfaces 
which emerges due to the FBC. Moreover, the use of 
PBC allows to obtain relative shifts of grains at Gl au- 
tomatically, as the coherent G2 grain boundary does not 
restrict the shift. 

Since Gl and G2 grain boundaries are not-equivalent, 
the following equation allows to calculate the specific 
excess energy related to the grain boundaries: 



■y = [EcB(n,a,c) - nE,uh{l,a,c)]IS , 



(1) 



where Ecsin, a, c) is the energy of the cell with GBs and 
n Ti atoms, Esub(^,a, c) is the energy of one atom in the 
commensurate Ti hep bulk supercell, and S = lac. The 
values of 7 for the considered grain boundaries before 
and after relaxation are given in Table [1] The energy 
decreases by more than two times after relaxation. 

3. 2. Considered positions of C, N and O atoms 

In order to study the interaction of impurities with 
Gl grain boundary, the same positions for C, N, and O 
atoms are used. We have considered five different po- 
sitions of interstitials within simulation cell: XI at Gl 
grain boundary (X = C, N, O) and X2-X5 at several 
distances from the Gl plane. The corresponding con- 
figurations of all atoms for mentioned positions are also 
named by X/ latter in the text. Locations of positions are 
shown within one relaxed simulation cell in Fig. |3] as 
the corresponding atomic configurations after relaxation 
are quite similar for all cases. The X2-X5 positions sit- 
uated in slightly deformed octahedral pores. The pore at 
XI position has more difficult topology than octahedral 
one. 

We made additional full relaxation to take into ac- 
count the possible influence of impurities on the grain 
boundary structure. Table |2] contains distance d from 
the Gl plane to the impurities after relaxation. The X5 
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Figure 1: (a) Structure of the supercell with S7 GBs after mathematical construction, (b) Structure of the supercell with Z7 GB after the full 
relaxation. The shift of grains is clearly visible. To illustrate the structure of the second GB, two simulation cells are shown. Atoms at GB and bulk 
regions are painted in different colours. The solid lines Gl and G2 denote grain boundary planes. 
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position at the center of two GBs is considered as the 
reference state of impurity in the bulk of titanium. 

One may note that distance between impurity and its 
periodic image in [0001] direction is quite small and 
some interaction is possible. We have checked that this 
interaction is only 10 % of the smallest calculated value 
of segregation energy. 



z 




Figure 2: The structure of G2 grain boundary. The marked atoms cor- 
respond to those in Fig. \T\ The additional atoms, obtained by periodic 
replication, complete the conventional cell of the tj-Ti phase in grain 
boundary. Here a and b corresponds to the lattice parameters of to-Ti. 
The z and y axes are coincide with [0001] and [1230] directions. 

3.3. Segregation energies and volumes of C, N and O 
atoms at a-Ti grain boundary 

The stability ofgrain boundaries depends on their for- 
mation energies [2]. It is well established that the reduc- 
tion of the GB specific energy due to the segregation of 
some impurities can significantly improve the thermal 
stability of nanostructure 1 10i l47ll . 

In order to determine possible decrease of GB energy 
due to the impurities, the segregation energies and vol- 
umes of C, N and O atoms at S7[0001](1230) CSL GB 




Figure 3: Atoms 1-5 indicate initial positions used for C, N and O 
within one relaxed simulation cell. 
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Distance from GB plane (A) 

Figure 4: The dependence of segregation energy E ,,,^ from the posi- 
tion of impurity. 

in a-Ti were calculated at K according to the follow- 
ing equations: 

E,eg = -Exi - Ex5 
Vseg = Vxi - VX5, 

where Exi and Vxi is the energy and volume of the sim- 
ulation cell with atom X at i position, Exs and Vx5 is 
the energy and volume of the same simulation cell with 
atom X occupying a grain interior site. The segregation 
energies and volumes are Usted in Table |2] The same 
trends can be derived for C, N and O. The segregation 
of impurities is unfavourable near Yl GB in all consid- 
ered positions. Moreover, the preference of positions 
is decreasing towards the Gl grain boundary plane (see 
Fig. m. The segregation energies for the same positions 
almost do not depend on the type of impurity, indicating 
that C, N and O interact with Yl GB identically in a-T\. 
The increase of segregation energy is accompanied by 
the increase of segregation volume. The increase of the 
volume of the cells is related to the expansion of the Gl 
grain boundary. In particular, the distances between Tl 
and T2 atoms is 4.23 A for C5, 4.32 A for C3 and 4.53 
A for CI configurations. 

The shifts of grains to each other are listed in Table |2] 
The values of shifts Ari in the [5410] direction are quite 
small (less than 5 % of titanium lattice constant) and do 
not deserve much attention. However, grains have no- 
ticeable additional shifts in [0001] direction at Gl plane 
due to the presence of impurities. The maximum shift 
is observed for X3 configuration, which is up to 30 % 
greater than Ar2 in pure GB. The dependence of shifts 
from the position of impurity is rather complex, how- 
ever in most cases the bigger value of shift corresponds 
to the bigger value of Eseg- 
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Table 2: Segregation energy E^eg (eV) and its separation into mechanical E,„ and chemical Ecj contributions. Geometrical data provided include 
lateral shifts of grains Ari in [5410] direction and Ar2 in [0001] direction (A), distance d (A) between impurity and Gl grain boundary plane, and 
segregation volume Vseg (h^) for all considered configurations. The C5, N5 and 05 configurations have the same parameters and denoted as X5. 
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To gain more physical insight on the similarity of 
segregation energies for different elements, we sepa- 
rated Eseg into mechanical and chemical contributions 
according to Geng et al. (l^ as follows. The mechan- 
ical contribution E„, was determined as the difference 
between the energies of the structures, where the impu- 
rity has been removed without subsequent relaxation of 
the host lattice. In other words, £,„ related to the Ti-Ti 
interaction induced by the impurity. The chemical con- 
tribution 

Ech ~ E^cg ~ E/jj (3) 

describes direct interaction between impurity and Ti 
atoms. The calculated values of £■„, and Ed, are listed 
in Table |2] For all configuration besides CI one can see 
the following regularity: (i) The chemical contribution 
dominates under £„, and correlates with Eseg', (ii) Val- 
ues of both contributions for different elements at the 
same positions are very close; (iii) The positive values 
of Ei„ for X2-X4 positions are related to the changes 
of Ti arrangement at GB and additional deformation of 
octahedral pores near the Gl grain boundary; (iv) The 
positive values of Ed, connected also with deformation 
of octahedral pores and depletion of charge density near 
the grain boundary. In general, such elements as C, N 
and O behave in a similar way. The segregation energy 
at C 1 position with negative chemical contribution is of 
particular interest and will be discussed below. 

To explain the positive values of segregation energies, 
we have studied electronic structure and distribution of 
charge density. The line profile of charge density be- 
tween TI and T2 atoms is shown in Fig. |5]only for the 
case of carbon. The further discussion will be made for 
X impurity because all three elements demonstrate the 
same behaviour Carbon, nitrogen and oxygen tends to 
form covalent spd-bonds with Ti atoms accompanied 
with the redistribution of charge density: accumulation 
between Ti-X atoms and depletion between neighbour- 
hood Ti-Ti atoms. For X2, X3 and in less manner for 
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Figure 5: Line profile of charge density between TI (OA) and T2 
(4.22 A) atoms for the case of carbon segregation at C5, C2, and C3 
positions. The profile for C4 is just below C5 and not shown. The 
CI profile is much higher than C5 profile due to the vicinity of carbon 
atom and also not shown. 



X4 such redistribution affects the GB region. Depletion 
of charge density between TI and T2 atoms results in 
a weakening of Ti bonds at Gl GB, increase of their 
length and increase of Gl grain boundary width. The 
increase of GB width and volume results in a higher seg- 
regation volumes and energies. The density between TI 
and T2 is the highest for XI position due to the vicin- 
ity of impurity, but the atomic size effect leads to the 
biggest width of GB and the highest segregation energy 
compared to other positions. 

Kwasniak et al. recently obtained from first princi- 
ples that C, N and O decrease formation energies of 
stacking faults in a-Ti. This means that segregation 
of these elements is energetically favourable at stacking 
faults. However, authors does not explain the physical 
reasons for such behaviour Stacking faults has smaller 
formation energies than high angle GBs. Trelewicz et 
al. f23\ within the framework of statistical thermody- 
namics showed that in binary polycrystalline system the 
segregation energy decreases (increases for choice of 



6 




Figure 6: Local coordination of interstitial sites at GB for CI (a), Nl 
(b) and 01 (c). There are eight Ti neighbours in the case of carbon 
and only six in the case of nitrogen and oxygen. 



sign of Eseg in ||25|]) with the reduction of GB forma- 
tion energy. Hence, the segregation of C, N and O can 
be possible on low-energy interfaces in a-Ti, but hardly 
feasible at high-angle grain boundaries with large for- 
mation energies. 

3.3.1. The features of grain boundary segregation 

The segregation energies for X2 and X3 configura- 
tions does not fit into the general picture. The impu- 
rity at X3 position, being slightly far (by 1 A) from GB 
plane than at X2 position, has higher segregation energy. 
The impurity at X3 position attracts T4 titanium atom 
(see Fig. [T]) reducing T3-T4 spacing. Grain bound- 
ary tends to conserve T3-T4 spacing causing additional 
shift of grains in [0001] direction (see Table|2]i. In turn, 
impurity at X2 position does not change arrangement of 
Ti atoms at GB and preserve the same shift as in ref- 
erence X5 position. In consequence, additional distur- 
bance of grain boundary structure by impurity at X3 po- 
sition results in higher segregation energy, than by that 
at X2 position. 

Nevertheless segregation of carbon at CI position (di- 
rectly at GB) is unfavourable, the negative chemical 
contribution in segregation energy is observed (see Ta- 
ble|2]). It means that the work needed to remove the car- 
bon while not permitting the Ti atoms to relax from CI 
position is larger by 0.29 eV than that from C5 position. 
The comparison of C 1 with N 1 and 1 shows that in the 
case of CI the Ti surrounding coordination is more com- 
plex and consists of eight Ti neighbours. The subtrac- 
tions of atoms from grain boundaries are show in Fig. 
|6] The local coordinations of Nl and Ol positions are 
distorted octahedral sites with six Ti neighbours. The 
restructuring of interstitial site in the case of carbon may 
be related to its higher valency and chemical capacity to 
form bonds comparing to N and O. However, due to the 
significant rearrangement of grain boundary structure, 
the loss in E„, overcomes chemical gain. 

To explain the reasons for the negative chemical con- 
tribution only in the case of carbon, we calculated site- 



projected partial density of states (PDOS). The PDOS 
of the interstitials (/5-orbitals) and surrounding Ti atoms 
(^f-orbitals) for XI and X5 configurations are shown in 
Fig. |7] The PDOS of titanium for XI configuration was 
obtained by averaging over the Ti atoms shown on Fig. 
|6] and for X5 configuration over the Ti atoms of octa- 
hedral coordination around the impurity. The most no- 
ticeable differences of PDOS between XI and X5 posi- 
tions are observed for carbon. For the majority of cases 
one may note the following tendency for PDOS which 
is consistent with the loss of energy for the case of im- 
purity at XI grain boundary site: (i) The shift of XI 
PDOS to the higher energies relative to X5; (ii) De- 
crease of density of states at several energies for XI 
case. However, there is qualitative difference for the 
p states of carbon comparing to oxygen and nitrogen. 
In the case of carbon the center of gravity of p states at 
XI site is shifted to lower energies (-0.067 eV), while 
centres of gravity of N and O p bands at XI are shifted 
to the higher energies (+0.124 eV and +0.113 eV re- 
spectively). This can be connected with the negative 
chemical contribution in the case of CI configuration. 

The smaller values of E„, for N and O is in consistent 
with the lower values of Vseg which means that N and 
O atoms require less space at the grain boundary than 
carbon. 

Summing up the results, it should be mentioned that 
the used method for Es^g calculation has one drawback 
related to the optimisation of GB structure. The draw- 
back is linked to the fact that it is hardly possible to 
obtain by first principles the shift of grains which corre- 
sponds to the global minimum of total energy. There is 
a possibility of relaxation of compared configurations to 
different local minima of energy. However, the obtained 
self-consistent picture of the interaction between impu- 
rity and GB allows us to claim that qualitative results 
will remain unchanged for the true ground states. More- 
over, the quantitatively correct results have little sense 
for the considered case, since experimentally measured 
segregation energies are always statistically averaged 
over a large number of different sites at general grain 
boundaries. 

3.3.2. Comparison with experimental data 

A few words should be said about agreement of our 
results with experimental observation of increased con- 
centration of C, N and O at GBs in Ti by Semenova 
et al. lfl6ll . Primarily, we do not claim that equilib- 
rium segregation of these elements is principally im- 
possible in a-Ti, since we considered only one type of 
grain boundary. However, there are also two interpreta- 
tions for the experimental data can be done as authors 
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Figure 7: Site-projected partial density of states of impurities and sur- 
rounding titanium atoms for XI and X5 positions. The Fenni level 
lies at the energy zero. 



do not give any. Semenova et al. provide atom probe 
tomography composition profiles of interstitial impuri- 
ties across the grain boundary after annealing at 623 K 
(see Fig. 4 in lfl6ll V On the one hand, the asymme- 
try and significant width of the observed concentration 
peaks are usually related to the non-equilibrium segre- 
gation due to the grain boundary migration ^ . On the 
other hand, the correlation between peaks for carbon 
and oxygen suggests that C-0 co-segregation had been 
observed. Hence, there is an intriguing possibility for 
grain boundary C-O co-segregation in ff-Ti without its 
independent segregation. 

Speaking about the possible mechanisms of thermal 
stabilisation in nanocrystalline pure titanium, the results 
of the present work are in favour of grain boundary 
pinning by bulk precipitates formed from solid solu- 
tion ifisli . 



4. Conclusions 

In summary, we have investigated interaction of C, 
N and O atoms with high-angle 27[0001](1230) grain 
boundary (GB) in a-Ti through the ab initio PAW poten- 
tial calculations. The specific simulation cell with two 
non-identical GBs have allowed to use the full periodic 
boundary conditions due to the coherency of the second 
GB. The first GB has the structure of unsymmetrical tilt 
boundary with lateral shifts, while the second GB repli- 
cates the structure of w-Ti phase with (000 1)„ || (1210)jj 
and [1230]„ || [0001]^^ orientation relationship. 

The segregation energies were calculated for differ- 
ent positions at different distances from GB plane. It 
was obtained that all considered elements has similar 
behaviour and prefer to be in octahedral site of bulk 
Ti. The maximum energy loss is observed for the po- 
sitions at GB interface. The analysis of geometry and 
electronic structure showed that the increase of energy 
is due to the redistribution of electronic density and re- 
duce of the bonding strength within the grain boundary 
area. 

Therefore, we believe that grain boundary pinning by 
small Ti-C particles is the main mechanism of thermal 
stabihty in nanostructured titanium OlSll . We are aware, 
however, that in a real system more complex processes 
such as grain boundary co-segregation and precipitation 
could occur, making further theoretical and experimen- 
tal investigations desirable. 
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